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Abstract 

We show how pointwi s e values of a function, f(x), can be accurately 
recovered either from its spectral or pseudospectral approximations, so that 
the accuracy solely depends on the local smoothness of f in the neighborhood 
of the point x. Most notably, given the equidistant function grid values, 
its intermediate point values are recovered within spectral accuracy, despite 
the possible presence of discontinuities scattered in the domain. (Recall 
that the usual spectral convergence rate decelerates otherwise to first order, 
throughout . ) 

To this end we employ a highly oscillatory smoothing kernel in contrast 
to the more standard positive unit-mass mollifers. 

In particular, post-processing of a stable Fourier method applied to 
hyperbolic equations with discontinuous data recovers the exact solution 
modulo a spectrally small error. Numerical examples are presented. 
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1 . 


INTRODUCTION 


Let f(x) be a bounded 27r-periodic function whose Fourier coeffi- 
cients are given by 

= 27 / f(y)e" lk ' y dy , -«. < k < » . ( 1 . 1 ) 

-TT 

It is well-known that whenever f is a smooth function its 
spectral approximation — consisting of the partial sums 

s N f W = f N (x) = l f (k)e ik * X , (1.2) 

|k|$N 

converges pointwise to f(x). A typical error estimate in this case 
asserts that for any x in the domain we have 

l f ( x ) " * c s l I f I I ( s ) *N S+1 , s > 1 . (1.3) 

Here and below, stands for (possibly different) generic constant 

bounds, and ||f||( s ) denotes the largest maximum norm of f and its 
first s derivatives, the maximum taken over the whole domain. 

We thus see that the decay rate of the truncation error on the 
left of (1.3) is restricted only by the degree of smoothness of the 
function f . In this sense, the spectral approximation is termed to 
be spectrally accurate . If, in particular, f is a C°°- function, the 
truncation error is rapidly decaying, faster than any fixed (= indepen- 
dent of N ) polynomial rate. Thus, the spectral approximation of C°°- 
functions enjoys the so-called infinite order of accuracy; this is in 
constrast to the usually slower convergence rate due to a fixed degree 
polynomial accuracy. 

Next, assume only the gridvalues f y = f(y^) are known, at the 
2N equidistant gridpoints y v = -it + vh , h = 2 tt/ 2N , v = 0,1, . . . ,2N-1. 
Invoking the trapezoidal rule, the (exact) Fourier coefficients in (1.1) 
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are approximated by discrete sums of these known gridvalues 


. 2N-1 -ik-y 

*(k) = 1 l f e v 

2N L v 


-N S k £ N . 


(1.4) 


The difference between the exact Fourier coefficients and their 
discrete approximation is also known to be spectrally small 


|f(k) - *(k)| $ C s ||f|| (s) N' S , s > 1 . (1.5) 

As a substitute to the (exact) Fourier coefficients appearing in the 
spectral approximation (1.2), f(k) , let us use their discrete counter- 
part, £(k) . The resulting new approximation is found to be exact at 

the gridpoints x = y . In other words, we arrive at the trigonometric 

(1) V 

interpolant . 

I f(x) = * (x) = f ^(k)e ik * X . (1.6) 

| k | *N 

The two type of errors committed in this case— the original trun- 
cation error in (1.3) padded with the aliasing errors in (1.5) — both 
are spectrally small. Hence, if f is smooth over the whole domain, 
then its pseudo-spectral approximation (1.6) is spectrally accurate 
even in between the gridpoints 

|f(x) - * N (x) I < c s Mf|l(- s ) N S+1 » S > 1 . (1.7) 

We also note that as in the spectral and pseudo-spectral cases 
(1.3) and (1.7) similar error decay is obtained with higher derivatives 
and in more space variables; the norm on the right-hand side of (1.3), 
(1.5), and (1.7) should be "raised" accordingly. Moreover, if the func- 
tion is in particular analytic , then the spectral accuracy is further 
improved to be exponential : let 2n > 0 be the width of analyticity 

strip with maximum modulus | |f|| ; then an error bound of the form 
C n l|f|| n e- Nri follows, e.g. [7], 

Unfortunately, the pointwise errors associated with the spectral or 
pseudospectral approximations, suffer from the limitation of being 
dependent on the smoothness of the function f over the whole domain 
(real or complex) and not just on its local behavior in the neighbor- 
hood of the point of interest. This dependence of the local conver- 
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gence rate on the global smoothness, which is reflected by (though not 
a consequence of) the error estimates (1.3) and (1.7), is indeed in- 
herent in both approximations. That is, the roughness of the function 
in one part of its domain decelerates the convergence rate in the 
smoother part of it. Most notably is the case of piecewise smooth 
functions: not only that Gibbs phenomenon is recorded jrt points of 

discontinuity, but in addition the spectral accuracy is lost at re- 
gions where the function is smooth. 

In this paper we show how pointwise values of the function f(x) 
can be recovered from the information contained in either its spectral 
or pseudo-spectral approximations, so that the accuracy solely depends 
on the local smoothness of f , that is, its smoothness in the neigh- 
borhood of the point of interest x . If, in particular, f is in- 
finitely smooth in that neighborhood, then the value f(x) is approx- 
imated within infinite order of accuracy. Most notably, we recover 
pointwise values within spectral accuracy, despite the possible pre- 
sence of discontinuities scattered in the domain. 

For such pointwise recovery, we should dismantle the above local- 
global coupling limitation, associated with the (pseudo-) spectral 
approximations. To this end, we employ a regularization kernel which 
is convoluted against the (pseudo-) spectral approximation in the 
usual fashion. Our regularization kernel consists of the product of 
two terms: first we introduce a cut-off function to localize the ker- 

nel in the spirit, advocated above; secondly, it is multiplied by the 
spectral approximation of the delta function (= Dirichlet kernel) , so 
that spectral accuracy is guaranteed. Convolution with the resulting 
kernel has then the effect of (locally) smoothing the spectral and 
pseudo-spectral approximations. 

The paper is organized as follows: in Section 2 we briefly dis- 

cuss those fundamentals of Fourier summation which will be later 
needed. Smoothing of the spectral approximation is described in 
Section 3. In Section 4, we similarly treat the pseudo-spectral 
approximation. It should be emphasized that the latter case directly 
involves only neighboring gridvalues , so that the construction of the 
pseudo-spectral approximation can be avoided altogether . In other 
words, (intermediate) pointvalues are recovered here via a locally 
supported yet spectrally accurate interpolation recipe. We remark 
that more general orthogonal families - other than the aforementioned 
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trigonometric one - can be used as well, to yield spectral smoothing: 
the notable examples of Legendre and Tchebyshev are brielfy sketched in 
Section 5. We conclude with numerical evidence which back up on theore- 
tical considerations. 

In [6] , Mock and Lax have shown how to recover within polynomial 
accuracy pointwise values of discontinuous solutions to linear hyper- 
bolic equations. They have employed a locally supported unit mass 
post-processing kernel with a finite number of vanishing higher moments. 

Our spectral smoothing is motivated by the Mock and Lax discussion — indeed, 
our regularization kernel based on the Legendre spectral approximation 
is intimately related to their kernel. Majda and McDonough and Osher 
[5] on the other hand, extending their previous study [4] with regard to 
the same problem, have employed a spectrally accurate smoothing proce- 
dure by operating directly in the Fourier space. Our smoothing in the 
real space rather than in the transformed one seems to offer more ro- 
bustness, resulting from the use of physical space localization; the 
latter is in fact the key element which enables us to apply our smooth- 
ing procedure to pseudo-spectral approximations. Moreover, it is also 
applicable* in conjunction with orthogonal families other than the trig- 
onometric one. 

This work has been motivated by the numerical studies of (pseudo-) 
spectral simulation of shock waves. However, in this paper we restrict 
our attention to the level of approximation only; applications to P.D.E. 
will be discussed elsewhere. 

2. PRELIMINARIES ON FOURIER SUMMATION 


Given a 27r-periodic function $ with Fourier coefficient 

A _ 1 / tt _ i k * v A 

$00 = (2tt) 'J_ n $( y) e 7 dy , its spectral approximation $ (x) is 


~ ,.J 2?/ <I, (y)e" lk ‘ y dye lk ’ X = / $(x-y)D h (y)dy ; (2.1a) 


kUp 


— 77 


here D (y) stands for the Dirichlet kernel 
P 


yy’ 


1 r ik*y _ _1_ sin (p+l/2)y 

2 * ikiip = 2u sin y/2 


(2.1b) 


normalized so that it has a unit mass 
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TT 

/ D (y) dy = 1 . (2.1c) 

-IT ” 

An (a' postriori) bound on the truncation error is given by 

I, ,1 Hk)e ik - X | * l |k| s |i(k)||k|- s 5 ||*|| l |k|- s s 
|k|>p |k|>p lSJ |k|>p 

? 11*11 ( S )P" S+1 s > 1 , (2.2) 

in agreement with (1.3), taking (4>,p) = (f,N) . Thus we have 

1/ 'Ky)Dp(x-y)dy-'Kx) I £ C s I I I (s) p_S+1 ’ s > 1 ; (2.3a) 

for later purpose, we quote here the special case x = 0, 

1/ < Ky)Dp(y)dy - $ (0) I $ c s l 1*1 I (s)P _S+1 » s > 1 . (2.3b) 

The above error bound is not the sharpest bound possible: let oj( ; ) 

denote the function* s modulus of continuity, then Kolmogorov’s result 

( 2 ) 

yields an asymptotri ally exact bound 

1/ *(y)D (x-y)dy-O(x) | ^ / u(^-;D S <5) sin^d? + o(p” S (o(hY 

-7T p ir p 0 P \ P / 

(2.4) 

Turning to the pseudo- spectral approximation, we have encountered the 
additional source of aliasing errors, due to discretization of the 
(exact) Fourier coefficients' integral. Invoking the aliasing rela- 
tions, e.g. [3], 

°o 

$(k) = l *(k+2jp) , -P * k $ p , (2.5) 

j=-~ 

(2p equidistant point’s interpolant is assumed). The aliasing errors 
do not exceed 

I *(k)-4>(k) | s l | *(k+2jp) | $ | |*|| f0 • l |k+2jp|" S ^ 
j^O [S) jVO 

$ C sl 1*1 I (s) p S ’ s * 1 * ( 2 -6) 
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in agreement with (1.5), taking (4,p) = (f,N). Hence, the aliasing 
error |$(k) - #(k) | , k = -p,...,p , adds up to a contribution 
similar to that of the truncation error, yielding in view of (2.5) 

l& p (x)-*(x)| .< | I $(k)e lk * X | + l |®(k) - &(k)| + i.|^(-p) + $(p)| $ 
y |k|>p |k|*p z 

* C S H*I I (s)P" S+1 » s > 1 , (2.7) 

in agreement with (1.7), taking (<K,p) = (f,N). It should be noted, 

(e.g. [3]), that there is no qualitative difference between the spectral 
and pseudo-spectral approximations. 

3. RECOVERING POINTWISE VALUES FROM THE SPECTRAL APPROXIMATION 


In this section, we show how to extract highly accurate approxima- 
tion to the point values of a discontinuous function from its first N 
Fourier coefficients in regions where the function is smooth. The basic 
idea is that these coefficients are moments of the functions and conse- 
quently, the integral of any smooth function against the spectral approxima- 
tion is highly accurate with that against the function itself. We 
therefore construct an auxiliary function such that when the spectral 
approximation is integrated against it, the desired original point value 
at a given point is recovered. 

To do that, let p(y) be a C S -function vanishing outside the 
interval (-it,tt) and normalized to take the value one at the origin 


p(y=o) = i . 


(3.1) 


We recall that the Dirichlet kernel in (2.1b) is given by 


V y) 5 2 ? 


|k| £p 

multiplying the two we obtain 


Y e ik * y = 1 sin (p+l/2)y . 

1 1 " 2t t sin (y/2) ’ 


(3.2) 


<l» 1 ,p (y) = p(y)D (y) 


(3.3) 


We now set as our regularization kernel 
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/ ,p (y) = e" V ,p (0 _1 y) = e‘ 1 p(e" 1 y)D p (e" 1 y) , (3.4) 

depending on a yet to be determined free parameter 6 , 0 $ 0 s 1 . 

A 

Given the spectral approximation, f N , we smooth its value via 
convolution with the above regularization kernel, computing 


f N * * 0,P OO = / f N (y')^ 0,P (x-y')dy' 


(3.5) 


In order to estimate the error, we decompose 


* * Q>V - f = (iff) * ^m’ P + * (* 6,P " + 


+ (f * tp 6,P - f) 


(3.6) 


The first term on the right vanishes in view of the orthogonality be- 
tween the N-degree (trigonometric) polynomial tp *” and the truncated 

A IN 

sum f., - f , 

N 


<f N - f) * *5 ,p - o 


(3. 7a) 


Thus we are left with two sources of error in this case: the truncation 

error in the second term 


t®’ p = (f N -f) * (^ e,p - i® >p ) 


(3.7b) 


and the regularization error in the third term, 


R 0,P = f * ij; 9 ,p - f . 


(3.7c) 


With regard to the truncation error T 9,p ; Young inequality implies 


ii(f N - f ) . (» e -p . 1 1 .< ni N -fii 

L 


(3.8a) 


and in the view of (1.3) we conclude that this term is spectrally small 


I T n’ P H 5 ll(f N -f) * W 8,P - (v' p )l I S 


(3.8b) 
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0 D 

Turning to the regularization error, R , we compute at a given 
fixed point x 

o o X+0TT , 

R (x) E f * ^ >P(x) ' fW = / f Cx') e P( 2 T L -)D p C^)dy'-fCx). 

(3.9a) 

Changing variables y = X ~ y ~ and making use of (2.1c), the regulariza- 
tion error is simplified into 


R (x) = f * ^ 9,P ( X ) " f 00 = / «> e,x (y)D (y)dy 

^ -7T * 

0 X 

where the auxiliary function $ ’ (y) is given by 


(3.9b) 


< J >0,X (y) = f(x-0y)p(y) - f(x) . 


(3.9c) 


In view of the normalization (3.1), <i> ,X (y) vanishes at y = 0 , and 

by appealing to the truncation error estimate quoted in (2.3b), we end 
up with 

I I - I / * e,x (y) D p (y)dyp c s ||* 0 ,x M (s) p “ s+1 . (3.9d) 


Added together, we have shown in (3.6) - (3.9) the following: 

Proposition 3.1 (Main Error Estimate) 

0 P 

Let \p ,F be the regularization kernel (3.4). Fix a point x in 
0 x 

the domain, and set $ ’ to be the auxiliary function in (3.9c). 

Then, the following error estimate holds 


f N * ii/ 0,p (x) - f(x) | <? c s Mf|MU e,p ll (s) N _s+1 


+ C S ,I$ 


e,xi 


-s+l 


(S) 1 


(3.10) 


The following two lemmas, whose technical proofs are postponed to 
the end of this section, provide us with the necessary explicit bounds 
on the two terms appearing on the right of (3.10). 

Lemma 3.2 

0 P 

Let ip denote the regularization kernel in (3.4). The follow- 
ing estimate holds 


(s) 


<: 6 


-s-1 



(1+P) 


s+l 


ll*"' p ll 


(3.11) 
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Lemma 3.3 

0 X 

Let $ * denote the auxiliary function in (3.9c). The following 
estimate holds 


||* 9 ’ X || * (l + 0) S .||p|| • Max | D k f (y) | . (3.12) 

(S) ( ' 5J |y-x|*0ir 

O^kss 

g 

Choosing p = N , 0 < 6 $ 1 , we conclude from (3.6), (3.8b), 
(3.9d) and the last two lemmas, the main result of this section, stating 


THEOREM 3.4 


Let f be a bounded 2iT-periodic function with a given N-degree 
spectral approximation f ^ . Setting the regularization kernel 


.O.N 6 . . 1 , -1 . sin (N B + l/2)y/6 

* (y) = 2^0 p(6 y) — s T^- ( 7 /2 9 ) ' * 


(3.13) 


then for any x in the domain , we have the pointwise error estimate 

g 

|f N * * e ’ N (x) - f (x) I .< C s ||p|| (s) .N B [N.0- S ||f N -f|| J.N -0 " 858 + 

L 


+ Max |D k f(y) | «N -Bs ] . 
| y-x | £0 tt 
O^k^s 


(3.14) 


Choosing 0=6=1 brings us back to the exactly same global error 

estimate we had in (1.3). Taking 6=1/2 on the other hand, the 

truncation and aliasing errors' contributions in (3.14) are balanced, 

and we are led to the following: 

Corollary 3.5 (Spectral Smoothing) 

2 s 

Let p(y) be a C - function , supported in [-iT,ir] and satisfy- 

ing (2.1). Then, for any x in the domain, the value f(x) can be 
“ — A 

recovered via the spectral smoothing of , which obeys the following 

error estimate 

|f N * * 6 ^(x) - f (x) I s C s -||p|| (2s) .[N.0' 2s .||f|| + 

+ Max | D k f (y) ]-N _s+1 , s > 1. (3.15) 

| y-x | £ 01 ; 

0<k^2s 
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In general, of course, the choices of the cut-off function, p , and 
the 6-exponent, 0 < 6 < 1 , provide us with a whole variety of admis- 
sible kernels, for which we have: 

Corollary 3.6 ( Infinite Order of Accuracy) 

oo 

Let p(y) be a C - function , supported in [ —it , tt ] and taking the 
value one at the origin . Assume the function f is C°° in the t- 
neighborhood of a point x in the domain . Then the spectral smoothing 


f * J; 0=£/7r ’ N6 (xl = — 
N V lXJ 2e 


X + G „ 

/ f N (y)p 

y=x-e 


(^) 


sin(N e ijiiiizl 

e 


sm 


”(x-y) 

2e 


dy 


0 < 6 < 1 , (3.16) 

recovers the function value f(x) within infinite order of accuracy . 
Remarks 

(i) Suppose f is known to be smooth in the asymmetric neighborhood 

of x , (x - e L , x + e^) , 0 < $ tt . Let p be a C -function 

supported in the interval [-6 ^e. , 6 ] inside of [ -tt , tt ] , such 

L K 

that p(y = 0) =1 . Then a nonsymmetric version of the above spectral 
smoothing reads 




f N * * 


e,rr, , i x+ r e R - r , ,x-y, sin ^ + 

(x) -2 5Fe _/ f N ( y )p( e 3 — 


x-e. 


L N' 


dy , (3.17) 


sin 


20 


recovering f(x) within spectral accuracy. The case c = e = e = it 0 

L R 

coincides with Corollary 3.6. 

(ii) The above error estimates concerning the spectral smoothing 
f^ * ip still enjoy the further flexibility in choosing different 
s-orders in (3.10). This provides us with even further richness so as 
to tune the different free parameters to yield accurate results. 

As promised, we conclude this section with the following: 

Proof of Lemma 3.2 

0 D 

With the regularization kernel ip in (3.4), Liebnitz's rule 


gives us 


k 6,P | l s? e s_1 | | p (y) D (y) 


(s) 


(s) 


£ 0 


-s-1 


l (•) 

j=0 2 


Cs-j) 11 p 1 1 (j) » 


(3.18) 


complemented by the maximum norm estimate 
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'V'tJJ 5 


_L y i k ij < L_ J * 1 

2 " |k|Jp 11 ' ’« +1 > p ’ 


(3.19) 


the desired result follows 

s+1 


* 


6>P| 


I * I e- s+1 
Vs) * 


.1 ftOUPII^/.e-'llpIl^a.p)- 1 . 

J “ (3.20) 


Proof of Lemma 3.3 

0 X 

Let $ ’ (y) = f(x-0y)p(y) - f(x) be the auxiliary function in 
(3.9c) with p (y) supported in [-ir,ir] . We observe that the only f- 

0 X 

values participating in the definition of 4> * are those from the 
9* it neighborhood of x, |y - x[ < 0-tt . Applying Liebnitz's rule 
restricted to that neighborhood, we find 


k0»x. 


« * ,1® 


l p l l m 0 S ~ J . M f x 1° f l 5 

| y-x | <©TT 

0<k<s- j 


$ (1+0) | |p | | * Max |D f (y) | 

^ J |y-x|<0TT 
0<k<:s 


(3.21) 


as asserted. 


4. RECOVERING POINTWISE VALUES FROM THE PSEUDOSPECTRAL APPROXIMATION 


In this section we treat the case where the discrete gridvalues 
f v = f(y v ) are given, so that a pseudo-spectral approximation 
collocating these gridvalues is uniquely determined; see (l.(>). The 
key observation here is that the integrand f^(y')^ ’^(x-y 1 ) in (3.5) 
is smooth over the whole domain, due to the kernel localization in the 
neighborhood of the point of interest, x . Hence, replacing the con- 
volution integral with an appropriate trapezoidal sum, only an addition- 
al spectrally small aliasing error is committed. Thus, in analogy with 
(3.5), we smooth the pseudo-spectral approximation via the convolution 


sum 


2tt 

2N 


2N-1 


l 

v=0 


^ (v )il> e, P(x-y ) = — 
1 y v J “ 2N 


2N-1 Q 

J f ip ’^(x-y ) . 
v=0 v 


(4.1) 


O 

Observe that since 4 1 is supported in the neighborhood of x , only 
those neighboring gridvalues are taking part in the pseudo-spectral 
smoothing. 


The computed error at a fixed point x , amounts to 
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w f v * 6 ' p ( x -y v ) - - 

v=0 


/ 2N-1 A \ / \ 

■ (§ l f/ >P Cx-y u ) - f * ♦ e,P w) * (f * * 8,P 00 - f(x)) . 

V=0 (4.2a) 

There are two sources of errors in this case: the aliasing error due 

to the use of the trapezoidal rule in the first difference 

2N-1 

A^ ,P = |£ 1 f y ,P (x-y ) - f * * e ’ P (x) , (4.2b) 

v=0 

and as before, the regularization error in the second difference 

f . (4.2c) 


R 0, P = f * /’P 


The aliasing error estimate in (1.5)^ and the regularization error 
estimate in (3.9d) yield: 

Proposition 4.1 (Main Error Estimate) 

Let *P be the regularization kernel (3.4). Fix a point x in 
the domain and denote 


X e,p,x (y) = f(y) •^ e,p (x-y) 

s 6 » x 


(4.3) 


Also , let 0 ’ be the auxiliary function in (3.9c). Then, the follow- 
ing error estimate holds 


„„ 2N- 1 

>2^ V f ii, 0 »P( 


2N V 

v=0 


C s llx°’ P ' I 'll (s) N- S * 


c s ll <t °’ x M C s)P " s+1 • 

(4.4) 


Q M y 

We observe that the newly introduced auxiliary function x (y) 

is supported in the 0 -Ti-neighborhood of x , where Liebnitz's rule 
yields 

n. e >p> x i 


CO * l cbl k 0,P | | fil - Max | D k f (y) | ; 

1 J j=0 J UJ | y-x | $6tt 

Osk$s- j 

invoking (3.11), the following bound is found 


(4.5a) 
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( e >P> x | 


(s) 


s 

I 

j=0 


(J)e- j+1 | 


(j) 


Ci+P) 


^ +1 « Max |D k f(y) 
|y-x|$07r 
O^k^s 


$ 0(l+p) | |p 1 1 , >• Max 

1 ’ | y-x | $0TT 
Oskss 


| D k f(y )|(ll|lE) 


s 


(4.5b) 


The last estimate on the aliasing part of the error, augmented with the 
previously derived estimate on the regularization error in Lemma 3.3, 
lead us to the main result of this section, stating: 


THEOREM 4.2 
Let f 


f 

v 


f(y v ) • 





(y) 


be a bounded 2-TT- periodic function with 
Setting the regularization kernel 


1 

2tt0 


p(e *y) 


sin (N^+l/2)y/0 
sin y/20 


given gridvalues 


(4.6) 


then for any x in the domain , we have the pointwise error estimate 


'2tt 2N y J 0,N P . 
I fj (x-y v ) 


2N 


f(x) 


v=0 


* C s I I P I Ifs) 


Max |D k f(y) |- [N" es+1 + 0" s .N B .N" (1_e:)s ] . (4.7) 

| y-x | ^0-rr 
0<:k<:s 


Taking 6 = 1/2 to balance the two error's contributions 1 , we find: 

Corollary 4,3 ( Pseudo-Spectral Smoothing) 

2 s 

Let p(y) be a C - function , supported in [-tt,tt] and satisfy- 
ing (2.1). Then for any x in the domain, the value f(x) can be 


recovered via the pseudo-spectral smoothing of the neighboring grid - 
values , f , which obeys the following error estimate 


I — 

1 2N 


n V 

v=0 


( x -y v ) 


f(x) 


* C s^ p ^(2s)‘, M ‘? x |D k f(y) | (1 + 0 _2 s )-N" s+1 , s > 1. (4.8) 

| y-x | S07T 
0£kj?2s 

In analogy with Corollary 3.6, we also have: 

Corollary 4.4 ( Infinite Order of Accuracy) 

oo 

Let_ p (y) be a C - function , supported in [-tt,tt] and taking the 
value one at the origin . Assume the function f is C°° in the e- 
neighborhood of a point x in the domain . Then the pseudo-spectral 
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smoothing 

2N-1 






2tt 

■2FT 


2N-1 0 / .3 v v /7r(x-y )\sin(N +-=■) 

y f * 6 = e/ ' i * n Cx . y ) = i- y f p f .l. V . ) — — s 

v=0 v J 2cl, y yx-e v ^ £ ' . "^-y, 

V sin 2e 


V ■ 


0 < 6 < 1 (4.9) 

recovers the function value f(x) within infinite order of accuracy . 

In closing this section, we would like to emphasize another, 
slightly more global variant of the pseudo-spectral smoothing, based on 
integral convolution of the pseudo-spectral interpolant against the 
regularization kernel 


y * ,0,n b _ Q* * /,0,n b T^e, n b \ 

*N * - f = V " V ) 


(y * * b,n b \ /. . .0,^ _\ 

V f N % - f * J + ( f * * - fj. (4. 


10) 


The first term is spectrally small due to the interpolation 
error associated with the smooth regularization kernel as argued in 
Section 2; by the exactness of the trapezoidal rule applied to (trigono- 
metric) polynomials of degree < 2N , we have 


y 'V© ?7 t 6 N® 

VXr =i I V 9 ’ (x-y v ) 

v=0 


(4.11) 


and consequently, the second difference is spectrally small as argued 
above in relation to the aliasing errors. Finally the third difference, 
is the spectrally small regularization error. 


5. CONCLUDING REMARKS 


The above arguments also apply to other orthogonal families. In 
conjunction with Legendre polynomials , we set as our regularization 
kernel 

if> 6,p (y) = 0 _1 p(e" 1 y)K p (e" 1 y) ; (5.1a) 

here p(y) is C S -function supported in the interval [-1,1] such that 
p(y=0) = 1 and K (y) stands for Legendre spectral approximation of 

the delta function 



-15- 


k Cy) 


= ? (k4: 


k=0 


2) p k (y) p k (°) 


(5.1b) 


normalized to have a unit mass 


/ K (y)dy = 1 . 

-1 p 

In view of the Christoffel-Durboux identity, we can rewrite 

, (r? - ki Vi C 0 )p pW - p p ( 0 >ViW 

P ’ 2 y 


(5.1c) 


(5. Id) 


The resulting spectral smoothing via the above Legendre-type regu- 
larization kernel was introduced in [1], and is intimately related to 
Mock and Lax [6] post processing: indeed, ^ 0,p serves as a locally 

supported kernel with vanishing higher moments and unit mass - modulo a 
negligible spectral error . 

Similarly, we can use Tchebyshev orthogonal expansion where 
K (y) in (5.1a) is replaced by 


K p (y) 



V y)T k C0) = 


d-y 2 ) 

jp TT 


y 


(5.2) 


We note that the (pseudo-) spectral smoothing done with the Tchebyshev 
kernel is not translated to the usual cut-off in the transformed space. 


6. NUMERICAL EXAMPLES 


In this section we demonstrate the efficacy of the smoothing pro- 
cedure outlined above. As a test function we have chosen the piecewise 
C°°- function 


f(x) 


• X 

sin 

. x 
(-sin j 


0 $ X £ IT 

TT £ X £ 2lT . 


( 6 . 1 ) 


As before, denote its spectral approximation by f^(x) , and let £^(x) 
be the pseudo-spectral approximation to f(x). It is evident from the 

A 

first column of Tables I and III that f^(y^) - the spectral approxima- 
tion sampled at y^ = vtt/N - do not approximate f(y^) within spectral 
accuracy. In fact, the error committed by f 128 (y ) is only half of 
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that committed by f 6 4 (y v ) > this is in accordance with a suitably 
sharpened error estimate of type (1.3) - consult e.g. (3.4). Regarding 
the pseudo-spectral approximation, i^(x) > it °f course collocates the 
exact values at the sampling gridpoints, i^(y y ) = f(y ) ; yet, in 
between these gridpoints, ^ N (y v+ i /2 = Cv + 1/2 ) tt/N) approximate 
f(y v+ l /2 ) within first order accuracy only, as shown in the first 
column of Tables II and IV. 

In order to construct our regularization kernel, we define the cut- 
off function P(5) = P a (5) to be 


exp 


P a (5) = 


? 2 -l 


U < 1 


otherwise 


( 6 . 2 ) 


oo I | 

namely P a (?) is a C -function whose support is the interval |£| < 1. 
Our regularization kernel is now of the form (see (3.4)) 


0 


,p (y) 


l 

2ir0 


p a ( 0 V) 


sin (p+l/2)y/0 . 
sin y/20 


(6.3) 


The post-processing procedure of the spectral approximation f. 


involves convoluting f^ against 


e,p 


namely 


2tt 


f(x) 


2tt0 


f f fy'jp f 2SlX) sin (p 1 l/2)(x-y)/9 . 
; r ^*.yj v 0 J sin (x-y)/20 


(6.4) 


where x is a fixed point of interest. (In practice we use the trape- 
zoidal rule to evaluate the right-hand side of (6.4) taking a large 
number of quaderture points.) 

The parameter 0 was chosen as 


0 = |x - ir| ; 


(6.5) 


this guarantees that ^ is so localized that it does not interact with 
regions of discontinuity. 

It should be noted, in this stage, that if 0 was so chosen to be 
the same for each x , (and not as in (6.5)), the formula (6.4) admits 
a simpler form; that is, if 
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/• P (y) = l o k e iky 

It-” 


then 


N 

l 

k=-N 


f(x) Vi f(k)a k e ikx 


( 6 . 6 ) 


(6.7) 


This procedure can be carried out efficiently in the Fourier space. 

Next, we turn to the post-processing for the pseudo-spectral 
approximation f N (x) which is simpler than (6.4). In fact, in this 

CaSG 2N-1 

f(x) I f (y v )^ e,P (x-y v ) . (6.8) 


Note that carrying out the smoothing procedure defined in (6.8) does not 
involve any extra evaluation of £(y) in points other than y^ , in 
contrast to spectral smoothing procedure in (6.4). As before, the para- 
meter 6 was chosen according to (6.5). We have yet to determine the 
parameters p and a . The parameter p must be equal to N^ for 
0 < B < 1 in view of (3.14), in order to assure infinite accuracy. (In 
our computations 8 % .8.) Finally we feel that a is problem depen- 
dent and we chose a = 10. We have not tuned the parameters to get op- 
timal results; further tuning may improve the quality of our filtering 
procedure. 

In Tables I, II, III, and IV we give the results of the smoothing 
procedure at several points in the domain. The pointwise values are 
now- recovered with high accuracy. The first column in each table indi- 
cates the points in which the procedure was performed. We limited our- 
selves to four points in the interval (0,tt) because of the symmetry 
of the function f(x) . 

The second column gives either the spectral approximation f (x) 
or the pseudo-spectral approximation ? N (x) , N = 128 in Table I and II 
and N = 64 in Tables III and IV. The third column gives the smoothed 
results, when filtered by (6.4) on (6.8), at the same points as in 
column I. 

The results indicate the dramatic improvement obtained by the 
smoothing procedure. Moreover, note that the error committed by £ 

r * . A 128 

(° r f 12 g) is better than the one committed by (or f fi4 ) only 

by a factor of 2 whereas after the post-processing the error 
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improves by a factor of 10 . 

Table I. 

Results of smoothing of the spectral approximation of f(x), N = 128 . 


TTV 

X = — o 

v 8 

If'v'-VVI 

l f - * 'H 

v equals 


at x = x 

V 

2 

3.2 (-3) 

5.8 (-10) 

3 

5.2 (-3) 

7.9 (-10) 

4 

7.8 (-3) 

6.3 (-10) 

5 

1.1 (-2) 

1.1 (-10) 

Table II. 



Same as Table I for the pseudo-spectral approximation ^(x). 

! 

TT * 

V)* " 8 t ' , * a ' 

v equals 

| f (x . ) - ? ,(x j ) | 
v+h N v+J* 

l f - * *1 

at x = x , 

v+h 

2 

5 (-3) 

7 (-10) 

3 

8.1 (-3) 

7.9 (-10) 

4 

1.2 (-2) 

6.4 (-10) 

5 

1.8 (-2) 

1.2 (-10) 

Table III. 



Results of smoothing of the spectral approximation of f(x), N = 

TTV 

X = "77 
V 8 

v equals 

l f (V> - Wl 

If - f N * *1 

at x = x 

V 

2 

6.4 (-3) 

4.8 (-6) 

3 

1 (-2) 

5.9 (-6) 

4 

1.5 (-2) 

7.7 (-6) 

5 

2.3 (-2) 

8.9 (-6) 
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Table IV. 

Same as Table III for the pseudospectral approximation, f ^(x) . 


Vi- 

v equals 

l f <v,’ - Vv’l 

If - * *1 

at x = x , 
v+*5 

2 

1 C-2) 

4.1 C-6) 

3 

1.6 (-2) 

6 (-6) 

4 

2.4 C-2) 

7.8 (-6) 

5 

3.6 (-2) 

8.9 (-6) 


7. ENDNOTES 

* The single and double priced summations indicate halving the first 
and the last terms, respectively. It is used in this case to com- 
pensate for the use of even number of gridpoints. 

2 

Referring to the convex case. 
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